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Abstract. We consider examples of the H-type groups with the natural horizontal distribu- 
lf~) • tion generated by the commutation relations of the group. In the contrast with the previous 

CN| ' studies we furnish the horizontal distribution with the Lorentzian metric, which is nondegen- 

erate metric of index 1 instead of a positive definite quadratic form. The causal character 

Ois defined. We study the reachable set by timelike future directed curves. The parametric 
equations of geodesies are obtained. 
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1. Introduction 



Sub-Riemannian manifolds and the geometry introduced by bracket generating distributions 
of smoothly varying £;-plains is widely studied, interesting subject, which has applications in 
control theory, quantum physics, C-R geometry, and other areas. The main difference of the 
sub-Riemannian manifold from a Riemannian one is the presence of a smooth subbundle of the 
tangent bundle, generating the entire tangent bundle by means of the commutators of vector 
"^j- ■ fields. The subbundle, which is often called horizontal, is equipped with a positively definite 

. metric that leads to the triple: manifold, horizontal subbundle, and Riemannian metric on the 

horizontal subbundle, which is called a sub-Riemannian manifold. The foundation of the sub- 
Riemannian geometry can be found in |17l fl8j fl9| I20| . The following question can be asked. 
What kind of geometrical features will have the mentioned triplet if we change the positively 
definite metric to an indefinite nondegenerate metric. It is natural to start with the Lorentzian 
metric of index 1. In this case the triplet: manifold, horizontal subbundle, and Lorentzian 
metric on the horizontal subbundle can be called sub- Lorentzian manifold. It was mentioned 
in [20] that it would be interesting to consider the sub-Lorentzian geometry, but, as it is known 
to the authors, there are only few works devoted to this subject [T], [TT1 [12] . In the present paper 
we study two examples of H-type groups furnished with the Lorentzian metric. The notion of 
H-type group was introduced in p3], see also a nice classification in [9]. The simplest example 
of a nonabelian H-type group equipped with the Lorentzian metric is the Heisenberg group with 
the Lorentzian metric and it has been considered in I12j . Our example is an extension of 
the Heisenberg group having 3-dimensional center related with the quaternion division algebra. 
The so-called quaternion H-type group with the positive definite metric was studied in [5], [5], [6] . 
We study the Heisenberg group and the quaternion H-type group endowed with the Lorentzian 
metric defining a causal character of the manifold under consideration. We give a description 
of the reachable set by timelike future directed curves. We find the parametric equations for 
geodesies on the Heisenberg and quaternion H-groups with the Lorentzian metric. Unlikely 
to the sub-Riemannian Heisenberg group, we get the uniqueness of geodesies starting from 
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the origin. Notice that by a geodesic we mean a projection of a solution of the corresponding 
Hamiltonian system onto the underlying manifold. 

2. Basic definitions and notations 

H-type homogeneous groups are simply connected 2-step Lie groups G whose algebras Q are 
graded and carry an inner product such that 

(i) Q is the orthogonal direct sum of the generating subspace V\ and the center V 2 : Q = 

V 1 ®V 2 , V 2 = [V 1 ,V 1 ], [V u V 2 ]=0, 
(ii) the homomorphisms Jz '■ V\ — ► V\, Z G V 2 , defined by 

(j z x,x') = (z,[x,x'}), x,x'ev u 

satisfy the equation j| = — \Z\ 2 I, Z G V 2 . 

Here (-, ■) is an inner product on Q, [•, •] is a commutator and I is the identity. 

We will consider the quaternion H-type group H with V\ associated with the space of 
quaternions and V 2 as a three-dimensional center. Remind that quaternion numbers are a 
noncommutative extension of complex numbers. A quaternion can be represented in a matrix 
way by q = aU + bX\ + cL 2 + dT 3 , where 

10 0" 
-10 
1' 
0-10 

"0 -1 o" 
1 
10 
0-100 

is the basis of quaternion numbers given by real matrices. As a vector space Q over the real 
numbers, the quaternions have dimension 4, whereas the complex numbers have dimension 2. 
H is a set M 4 x R 3 with the multiplication law defined in the following way 

L qW) = L {x,z){x' ,z') = (X,Z!,Z 2 ,Z 3 ) O ( X ' , z[, Z 2 , z' 3 ) 

= (x + x',zi + z[ + ^(lix,x'), z 2 + z' 2 + ^(2 2 x,x), z 3 + z' 3 + i(J 3 x,x')), 

for q = (x,z) and q' = (x',z'), where x,x' G M 4 and (Xix,x'), (T 2 x,x'), (l 3 x,x') are the usual 
inner products of the vectors T\x, T 2 x, l 3 x belonging to R 4 by x' G M 4 . The unit element e 
of the group H is e = (0,0) and the inverse element to (x,z) is (— x, — z). The multiplication 
"o" defines the left translation L q of q' by the element q = (x,z) G H on the group H. The 
topological dimension of the group H is 7, the Hausdorff dimension is 10. The matrices 1\, 1 2 
and Z3 define the homomorphisms Jz- 

Q is a Lie algebra of the group H which can be associated with the set of all left invariant 
vector fields of the tangent space TH e at the unity e. T^H C TH is a subbundle of the tangent 
bundle spanned by the left invariant vector fields X a (x, z) with X a (0, 0) = tJ-, a = 1, ... ,4. 

In coordinates of the standard Euclidean basis a = 1, . . . , 4, tJ-, (3 = 1, 2, 3, these vector 
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10 
10 
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fields are expressed as 

X - 9 + 1 

9 1 ( 

9 1 ( 
9 1 ( 

We call T/jH the horizontal bundle and any vector field V £ T^H is called the horizontal vector 
field. The left invariant vector fields Zp = gf^-, = 1,2,3 form a basis of the complement 
to T^H. The commutation relations are given by 

[Xi,X 2 ] = —Zi, [Xi,Xs]=Zs, [Xi,X^\ = Z2, 
[X 2 ,Xs] = Z 2 , [X2, A4] = —Z3, [XfaXj] = —Z\. 

We see that X a , a = 1, . . . , 4, and their commutators Zp, (3 = 1, 2, 3, span the entire tangent 
bundle TH. Therefore, the horizontal bundle T^H is bracket generating of step 2, see the 
definition of bracket generating subbundle in [SJ [T3] . 

A basis of one- forms dual to X\, . . . , X4, Z\,Z 2 , Z3 is given by dx\, . . . , dx^, vi,v 2 ,V3 with 

v\ = dz\ — - (+x 2 dx\ — x\dx 2 + xt±dx3 — £3^x4), 

v 2 = dz 2 — -(—x^dxi — X3dx 2 + x 2 dx3 + xidx±), 

V3 = dz3 — -(— x 3 dx\ + x^dx 2 + x\dx3 — x 2 dx^). 
We use the normal coordinates (x, z) = (xi, . . . , Zi, z 2 , Z3) for the elements 

/ 4 3 \ 

exp ^2 x aX a + zg-Z^ G H. 
\a=l p=l J 

An absolutely continuous curve c(s) : [0, 1] — > H is called the horizontal curve if its tangent 
vector c(s) belongs to V\ at each point c(s). Any two points in H can be connected by 
piecewise smooth horizontal curve because of bracket generating property of IX H. Unless 
otherwise stated all vectors and curves are supposed to be horizontal. 

The geometry of H-type groups was studied, for instance, in [U [21 El El El El ESI US]- In all 
these cases the horizontal bundle T^H was endowed with a positively definite metric. We will 
consider the same type of the horizontal bundle but equipped with the Lorentzian metric. We 
call such a manifold sub-Lorentzian, analogously to the sub-Riemannian case, and we formulate 
it in the following definition. 

Definition 1. Let M be a smooth n-dimensional manifold, let IX M be a smooth fc-dimensional, 
k < n, bracket generating subbundle on TM, and letQr A Af(', •) be a smooth Lorentzian metric 
on ThM. Then the triple (M,ThM,Qx h M '(*>•)) is called the sub-Lorentzian manifold. 

We will call Qt h m{', •) the sub-Lorentzian metric and skipping the subscript, we write Q(-, •). 
We define a smooth Lorentzian metric on T^H by 

(2.1) Q(X 1 ,X 1 ) = -1, Q(X a ,X a ) = l, a = 2,3,4, Q(X h X-) = if i + j. 
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Given a sub-Lorentzian metric Q x at x G H, we can define a linear mapping g x : T*H X —* 
TH X as follows: for £ G T*H X , the linear mapping Y — > £(Y) for 1" G TftH x is represented 
uniquely by the identity 

(2.2) Qx(y,fe(0) = Z(Y) for all F G T h H x . 

The properties of Q x implies that g x varies smoothly in x, is symmetric, and nondegenerate. 
If Q is of index 1, then index of g is also 1. 

Conversely, given a symmetric nondegenerate linear operator g x : T*H X — > TH X with the 
image ThH x , there is a unique nondegenerate quadratic form Q x satisfying (|2.2p . 

Example. Let us consider the sub-Lorentzian metric Q x and co-metric g x on the quaternion 
H-type group: 

/-l 0\ 
10 
10 
\ 1/ 



{Qx} 



a/3 



a/3 
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-1 













-\x 2 




|x 4 




|x 3 
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2 X4 


-5X3 




2 , 2 , 2\ 

- a; 2 + x 3 + :r 4 ) 




5X2X4 




5X2X3 








-§X 3 




§xi 




5X2X4 




1 2 | 2 2\ 
-f- X2 + X :i — X4 ) 




— 5X3X4 
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\X 3 


±x 4 


|xi 


-|x 2 




\x2X3 




— 5X3X4 




1 2 2 

+ x 2 -x 3 


+ xl)J 



The latter matrix is symmetric, thus, it can be diagonalised to 
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Therefore, the index of the co-metric g^ is 1 similarly to the index of the initial quadratic 
form Q. 

Let us define the causal character of H. Fix a point p G H. Denote by Sl p the set of all 
horizontal curves c(s) : [0, 1] — ► H starting from c(0) = p. A horizontal vector v G T^Hp is 
called timelike if Q(v,v) < 0, spacelike if Q(v,v) > or v = 0, null if Q(v,v) = and v 7^ 0, 
nonspacelike if Q(y, v) ^ 0. A horizontal curve is called timelike if its tangent is timelike; 
similarly for spacelike, null and nonspacelike curves. 

By the time-orientation of a sub-Lorentzian manifold (H, T^H, Q) we mean a continuous 
horizontal timelike vector field on H. The time orientation of (H, T^H, Q) is given by the 
vector X\ according to (|2.ip . Then a nonspacelike v G T^SLp is called future directed if 
Q(v,Xi(p)) < 0, and it is called past directed if Q(v,Xi(p)) > 0. Throughout this paper f.d. 
stands for "future directed", t. for "timelike", and nspc. for "nonspacelike". 

For an open set U and fixed p G U, we define two reachable sets: I + (p, U) (resp. J + (p, U)) 
is the set of all points q G U that can be reached from p along a t.f.d. (resp. nspc. f.d.) curve 
contained in U. In the Lorentzian geometry I + (p,U) is called the chronological future of p 
(with respect to U); similarly, J + (p, U) is called the causal future of p (with respect to U). 
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We define a length of a nspc. curve c: [a, /3] — > H in the following way 



L(c) 



\Q{c,c)\idt. 



Let ip: U — > R be a smooth function on an open set £/ in H. The horizontal gradient V^</? 
of <^ is a smooth horizontal vector field on U such that for each p G [/ and w G T^Hp we have 
(d v ip)(p) = Q(v,V Hf(p))- Locally it can be written in the following way 

4 

V h <p = -{X W )X 1 + ^{Xj^Xj, 
i=2 

where Xi, . . . , X^ is an orthonormal frame of T^H defined on [7 with X\ timelike. 
The following proposition can be found in 

Proposition 1. For any normal neighbourhood U of a point p 

(i) J + (p, U) is a closed subset in U ; 

(ii) cl(I + (p, U)) = J + (p, U), where cl stands for the closure with respect to U . 

3. Reachable sets for Quaternion H-Type Groups 

In this section we describe the reachable set for quaternion H-type group. We start from 

the case p = and / + (0) = / + (0,H). Thereto let us consider the family of functions rj a = 

4 

—x\ + x\ + x\ + x\ + a(|zi| + \z-i\ + | ^3 1 ) , \a\ ^ —= that represents the homogeneous norm 

V3 

on H. We want to describe the reachable set in terms of values of the function r\ a . We present 
the calculations for the positive values of z±, Z2, z 3 , other cases considered analogously. The 
horizontal gradient of n a is 



+ 



2xi 



\oi{x 2 



X 3 - X4 



Xl + 

x 3 + 



2x 2 + 2 a ( — Xl ~ X 3 + x i) 



Xo 



2x 4 + -^ a ( Xl ~ X2 



X3) 



X A . 



Consider 770 and the set 

(3.1) r = {m < 0, X! > 0}. 

Theorem 1. Let / + (0) denote the reachable set for H and Tq be as in (|3.ip . Then 

/ + (o) c r 

Proof. Since QCVhVo, X\) = —2x\ and (^(V^r/o, V^r/o) = 4ryo, we conclude that 770 is t.f.d. in 
To- We claim that 

(3.2) J+(0)cr = {7 ?0 ^0,xi^0}, 

where we set J + (0) = J + (0, H). Let p G J + (0). It means that there exists c(t) nspc. f. d. 
c: [0, T] — > M7 , c(0) = 0, c(T) = p. Since c is future directed it follows that 

Q(c(t),X 1 (c(t))) = -xi(c(t)) < 0. 

Hence the function xi(c(t)) increases from xi(0) = along c(t), and therefore, x\{c{T)) = 
x\{p) > 0. The nonspacelikeness of c implies that 

Q(c, c) = —x\ + x\ + x\ + x\ ^ 0. 
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We need to show that x\(T) ^ x\ (T) + x| (T) + x\ (T) , taking into account that x\(0) 
£2(0) = £3(0) = 2:4(0) = 0. It will be shown if we prove the inequality 

,3.3) + m + m > f t (^4^4) = Z2 ' Xi ■ 

V 2 "r" 3 ' 4 

since in this case 

d 



xi(t) > — [\]xl + xl + xlj and x\ (0) = (y as} + x| + xj) (0) = 



=> xi(T)>(yx| + x| + x2)(T). 

Consider the case X2 • £2 + ^3 • £3 + £4 • £4 ^ 0, because otherwise the inequality (|3.3p is obvious. 
Squaring both parts of (|3.3[) and simplifying, we get that (|3.3p holds if and only if 

(±2X3 - X3X2) 2 + (X2^4 - X 4 X 2 ) 2 + (^3^4 - X4X3) 2 ^ 0. 



Statement (13. 2h is proved. 

Now we look at the behaviour of any nspc.f.d. curve which projects onto the set 



x\ = \lx\ + x\ + x\, X\ > 



Due to the chain of inequalities 

d 



il = lit IV X 2 + X 3 + x 4j < + *3 + *4 ^ 

such curves must be null curves. In addition, since VhVo ls nspc.f.d. on To, we see that for 
any t.f.d. curve c: [0,T] — > M7 , c(0) = 0, the derivative along the curve c(t) 

^%(c(t)) = Q(V h rio(c(t)),c(t)) < 0, 

and the function t — ► ??o( c (^)) is decreasing a. e. It follows that I + (0) CTq. □ 
Corollary 1. Following the notations of the previous theorem we have 

r n {zi = o,z 2 = o, z 3 = 0} = j+(o) n = o, z 2 = o, z 3 = o}. 

Proof. Consider the straight lines 

' x\ = t cosh ip, 
X2 = t sinh ip sin ip cos 0, 
3:3 = t sinh y> sin ^ sin 9, 
^4 = t sinh 99 cos ^, i > 0, 



(3.4) 



where ip E [0, 27r], G [0, 7r] , y> G [— 00, +00] are constants. Any of these straight lines is a t.f.d. 
curve. They fill up the interior of To n {z\ = 0, Z2 = 0, Z3 = 0}. So, we conclude that 

r n { Zl = 0, z 2 = 0, z 3 = 0} c /+(o) n = o, z 2 = o, z 3 = o}. 

The inverse inclusion follows from Theorem [TJ □ 

Let us calculate sets, where ^7 hJ]a is nspc.f.d. We have 

a \/3a 
Q(V h T] a , Xi) = 2(-xi + -(x 2 - x 3 - x 4 )) = 2(-xi + ("y • n")), 
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Q(^hVa,^hVa) 



■1 + 16")^ 



x 2 



T 2 
x 3 



4) 



+ 



a 
~2 



((x 2 + x 3 ) 2 + (x 2 + x 4 ) 2 + (x 3 - x 4 ) 2 ) 
2 3a 2 



= 4((-l + — )(x 2 -M 2 ) + — I 

where we used the notation ~y = (x 2 ,xs,X4), it = -4^(1,— 1,-1) and " x " is the vector 

product in IR 3 and | ■ | is the length of the vector with respect to the inner product in R 3 . We 
see that VhVa is f.d for \a\ < in IV Indeed, if we write the condition of future directness 
as 



a 

-xi + -(x 2 



x 3 ~ x A ) = {N, x) < 0, 
then it holds if and only if 



N = (1,^,—,--), ~x = {xi,x 2 ,x 3 ,x 4 ), 



(N,N)<0, <0. 
The latter two conditions satisfy in To for \a\ < -j=. 
We introduce the notations 

r a = {v a < 0, x x > 0} 

and 

A a = {(-l + ^r)(4 - 1 ~y\ 2 ) + ^H~y x n\ 2 < o, x x > o}. 

If a = 0, then Aq = Tq, otherwise A a C Tq. If a = -4= then ij_ = {0} and Aj_ = {x 2 

— X 3 = — X4, Xi > 0}. 

Remark 1. The set vl a has the following geometric sense. Let 

/-l 0\ 

Q = {Qij} = 



10 
10 
\ 1/ 



let 



V = {V kl } = 



1 ( rp 2 ,-y>2 I I ,vi 2 \ 

4 x 2 ' x 3 "I" x 4/ 



\x 2 X4 



9X2X4 
I, 



ryX 2 X 3 
-5X3X4 



1 Z' •> ' 2 _[_ , . , 2 ,-y»2 1^ ^y.2 \ 

4^x 1 i-x 2 x 3 -t-x 4 j / 



4 (^l + x 2 + x 3 ^l) 
^X 2 X 3 -ryX 3 X 4 

be quadratic forms and let 1? = (xi, x 2 , x 3 , X4), 2 = (zi,z 2 ,zs) be vectors. Then A Q can be 
written as 

A a = {Q( x , ~x) + V(z , z) < 0, x\ > 0}, with zi = z 2 = z 3 = a. 
We notice that both of the quadratic forms appeared in the co- metric g a @. 
Theorem 2. We have 

i + (0)nA Q cr a . 

Proof. Let p G 7 + (0) n A a . Then there is a timelike f.d. curve c : (0, T) -> r n A*, c(T) = p. 
Since the vector V^t/q is timelike and f.d. in the same set we conclude that the derivative 
4fTja(c(t)) = Q(S7 hVonc) is negative and the function rj a (c(t)) decrease along c(t), that yields 
Va{p) < jointly with x\{p) > 0. □ 

The complete description of the reachable set / + (0) is complicated, therefore we study its 
subsets. 
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Proposition 2. Let B be one of the following sets {x% = x^ = 0}, {x2 = £3 = 0}, {X2 = £4 = 
0}. Then 

/ + (0)n5 = r 4 nB. 

Proof. If X3 = X4 = 0, then for any horizontal curve the horizontality conditions 

1, 

Zl = -{+X 2 X 1 - X X X 2 + X4X3 - X3X4), 

(3.5) z 2 = -(—X4X1 - x 3 x 2 + X2X3 + X1X4), 

^3 = -(-X3X1 + X 4 X 2 + X1X3 - X 2 X 4 ) 

implies that z\ = |(+X2Xi — X1X2) and Z2 = zz = 0. Since 21 (0) = z 2 (0) = 23(0) = 0, we 
conclude that this case is reduced to the Heisenberg sub-Lorentzian manifold and we can apply 
results of [IT] . Other cases are obtained analogously. □ 



4. Hamiltonian and geodesics in Heisenberg group with Lorentzian metric 

In this section we study geodesics on the one dimensional Heisenberg group equipped with 
the Lorentzian metric. We remind that the Heisenberg group H 1 is the space M 3 furnished 
with the non-commutative law of multiplication 

(x,y,z)(x',y',z) = (x + x' , y + y' , z + z + ^(yx' - xy)) . 

The two dimensional horizontal bundle T^H 1 is given as a span of left invariant vector fields 

did did d 

X = -^- + -yTr, Y = -^---x-^-, [X,Y\ = Z = i-. 
dx 2 dz dy 2 dz ' dz 

We suppose that the Lorentzian metric Q is defined on T^M 1 by 

Q(X,X) = -1, Q(Y,Y) = 1, Q(X,Y)=0. 

The time orientation is given by the horizontal vector field X. Thus the triple (M 3 , T^M 1 , Q) 
is called the Heisenberg group with the Lorentzian metric, and to differ it from the classical 
case we use the notation M\. The reachable set for the Heisenberg group with the Lorentzian 
metric was studied in [JTJ [12] . 

To study geodesics we apply the Hamiltonian method. The Hamiltonian is defined as a 
symbol of the corresponding wave equation formed by the left invariant horizontal vector 
fields. The definition of geodesics in sub-Riemannian geometry differs from the definition of 
geodesics in the Riemannian geometry. Geodesics in the sub-Riemannian geometry (and we 
shall use the same definition) are projections of solutions of the corresponding Hamiltonian 
system onto the underlying manifold. The interesting feature of the sub-Riemannian geometry 
that even locally there is no uniqueness of geodesics. The study of number of geodesics on 
Heisenberg group with positively definite metric can be found, for instance, in [21 [10]. We 
apply the method of Hamiltonian mechanics in order to calculate geodesics of different causal 
characters and to study the uniqueness problem. 
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The governed operator is an analogue of the sub-Laplacian which has the form 
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dx 2^ dz 
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4 y dz 2 
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' dx dz 



dx 2 dy 2 J 4 



2 i 2 
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d 2 



d 1 d 

dy 2 dz 
d 2 
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dz 

d_d_ 

^ dx dz 



d_d_ 

dy dz 
d d 



x- 



dz 2 dx dz dy dz 
Then the associated Hamiltonian function H (£, rj, 9, x, y, z) becomes 



H= 1 -(-e+r 1 2 )- 1 - 



(- x 2 + y 2 )6 2 + l -(-yt-xr ] )e, 



where we use the notations £ 



(4.1) 



~§x ' ^ = "By ' ^ = ~5i ' corresponding Hamiltonian system is 
dH . yO 



di - 

dH 

dr] 

dH __ 

~d~e ~ 

dH 

dx 



-e ■ 

_l(_ x 2 +y2)0+ l ( _^ 

= -\x6 2 + ^rt6, 



xrj), 



dH 

dz 



0. 



The initial data is x(0) = y(0) = z(0) = 0, £(0) = £ , r?(0) = r? , (9(0) = 6. We are interested 
in finding the projection of the solution of (|4.1|) onto the (x, y, z)-space. We reduce the sys- 
tem (|4.ip to the system containing only (x, y, z) coordinates withthe initial data zero. If we 
express £ and r/ from the first two equations and substitute them in the expressions for £ and 
77, we obtain £ = ^y0, fj = — ^x9. Taking into account that 9 is constant, we differentiate the 
first two equations and replace £ and r] there. We get 



or 



y 



-9 



We are looking for the solution x = x(t), y = y(t), t £ [—00, +00], satisfying x(0) = 0, 
2/(0) = and x(0) = xq = — £0, y{0) = yo = Vo- Eigenvalues and eigenvectors for the matrix 

M :— ( 5 ^ I are A = ±\9\ and ( \ V I 1 ) respectively Therefore, 



9 



1 




ciel 9 !* + c 2 e"i c 
— cie'^l* + c 2 e~ 



xq = ci + c 2 , 



-ci + c 2 . 
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It follows that 



(4.2) 



Finally, we obtain 




(4.3) 



x(t) 

y(t) 



xq cosh(|#|t) — yo sinh(|0|i), 
-x sinh(|0|t) + y cosh(|0|i). 



^ s i n h(|0|i)-^(cosh(|0|*)-l), 
-^(cosh(|^)-l) + ^smli(|^). 



Lemma 1. // a geodesic is timelike future directed (past directed) at t = 0, then it remains 
timelike and future directed (past directed) for all t G [0, ±oo]. 



Proof. Since geodesies satisfy the Hamiltonian system (|4. 1 j) . the speed is preserved along 
geodesies and this implies that the causality character does not change. This also follows 
from (|4.2p because — x 2 (t) + y 2 (t) = —x\ + y^. We remind that a geodesic is future directed 
if x(t) > 0. Let us show that if we start from xq > and — xq < yo < xq, then the ge- 
odesic remains future directed for all t G (0,±oo). If < yo < xq, then we conclude that 
x = xocosh(|#|i) — yosinh(|0|t) > for t > because of the inequality sinh(|0|i) < cosh(|0|i). 
If — xq < yo < then we take t £ (0, — oo) and get future directed geodesies because of 
— sinh(|^|f ) < cosh(|#|t). We conclude that for xq > 0, and |yo| < xq we obtain future directed 
geodesies for t G (0, ±oo). □ 



Using the condition of horizontality, we get 



z = \iv x - x v) 



1 

W\ 



-i;2 + y o 2 )(l-cosh(|0|t)). 



Notice that the derivative i is positive for timelike geodesies, for lightlike curves we have z = 0, 
and spacelike geodesies satisfy i < for any values of t. Integrating, we get 

|2 



(4.4) 



z(t) 



\VQ\ 



26 2 



(|0|t-sinh(|0|t)), 



Pol 



■ 2 , -2 

-^0 + 2/0- 



The value of z-coordinate for timelike future directed geodesies starting from the origin is 
positive for positive value of t and negative for negative value of t, lightlike curves has vanishing 
^-coordinate and spacelike curves admits the negative values for t > and positive for t < 0. 
Further, introducing the notation \\h\\ 2 = —x 2 + y 2 we get 



(4.5) 



2 i 2 

-x +y 



4[H 

e 2 



■sinh 2 (— ). 



Lightlike geodesies are projected into the set — x 2 +y 2 = and spacelike geodesies are projected 



If a geodesic is timelike, then its projection onto (x, y)-plane lies inside the domain — x 2 +y 2 < 0. 
Lightlike geodesies are projected into the set - 
into the domain — x 2 + y 2 > 0, see Figure [TJ 

Lemma 2. Let j(t) = (x(t), y(t), z(t)} be a geodesic. Then 



-x 2 (t) + y 2 (t) + i\z(t)\ 



< if j(t) is timelike, 
= if 7(t) is lightlike, 
> if j(t) is spacelike. 
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Figure 1. The graph of timelike geodesic passing through the origin 
Proof. We use the notation ||/i(t)|| 2 = — x 2 (t) + y 2 (t). If z > 0, then the value 
\\h(t)\\ 2 + A\z{t)\ = \\h(t)\\ 2 + 4z(t) = ^A-i + \e\t + e-W) 

is negative for timelike curves and positive for spacelike curves, because -1 + \9\t + > 
for all t 0. The value of \\h\\ 2 + 4z vanishes for lightlike curves. If z < 0, then the expression 

2IMI 2 , , 



+ 4|z(t)| 



4z(t) 



9 2 



is still negative for timelike curves and is positive for spacelike geodesies, because —1 — \9\t + 
e l l* > o for all t ^ 0. The value \\h\\ 2 — 4z vanishes for lightlike curves. □ 



Lemma 3. The equation 
(4.6) 



4z 



-x 2 + y 2 



sinh (r) 



coth(r) 



has a unique solution r for given (x,y,z) if —1 < 2 < 1> an d has no solution otherwise 



Proof. The function /x(r) 



sinh (t) 



coth(r) is strictly decreasing in the interval (—00, +00) 



from 1 to —1, see Figure [2j It proves the lemma. 




Figure 2. The graph of function /z(t) 
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□ 

We would like to draw the reader's attention to the function \x. It is intriguing analogue of 
the function used by Gaveau [10] to study geodesies and geometry in general on the classical 
Heisenberg group. The classical counterpart Jl of the function [i has the form Ji(t) = ^ ^ — 
cotr, and the classical Heisenberg analogue of the equation (|4.6|) . 

A\z\ T 

7) 77 — o COt T 

x z + y z sin r 

has more then 1 solution. 

We state the following theorems describing the reachable set by geodesies starting from the 
origin. 



4|«| 

Then there is a unique timelike future directed (past directed) geodesic, joining O = (0, 0, 0) 
with the point A. Let 9 be a solution of the equation 

Az _ \6\/2 

-x 2 + y 2 ~ sinh 2 (|#|/2) 

Then the equations of timelike future directed geodesic 7 : [0, 1] — > are 



Theorem 3. Let A = (x,y,z) be a point such that x > (x < 0),—x 2 + y 2 < 0, < 1. 

n y 



x(t) = sinh 2 (f t)^(coth(f t)coth(M) - 1) +y(coth(JJt) - coth(f 
y(t) = sinh 2 (f t) (y( coth(f t) coth(i|) - l) + x( coth(f t) - coth(f )) j . 



(4.8) 



(4 9) z(t) = Z MZ^M 

^ * W |0|-sinh(|0|) 

The square of the length of these geodesies is 

|||*(1)|| 2 | + 4|*| 



I' 



2(\\6\ - sinh(|6»|)| + 2sinh 2 (|6>|/2)) ' 



Proof. Let us suppose that geodesies joining O = (0, 0, 0) with A are parametrized on the 
interval [0, 1]. We fix the solution of the equation (j4.7j) for a given point A = (x, y, z), x > 0, 

II II 2 

— 1 < x ^lyi < 1- Put t = 1 in (|4,4p and find = rarr^gngj ■ Then we substitute the last 

expression in (|4.4p and get (|4.9p . 

To obtain (|4.8p we do essentially the same. Set t = 1 in the expressions for x(t) and y(t) to 
find xoi 2/0- We get 

^0 = y(xcoth(y) +y), y = y(ycoth(y) +x). 

Substituting these expressions in (|4.3p . we get (|4.8p . 

Let us calculate the length of timelike geodesies. Expressing ||t>o|| 2 from (|4.5p . we get 

||2 _ \\h(l)\\ 2 9 2 
" °" 4sinh 2 (f)' 

Then the length of geodesies is 



(4.10) ^r^a^vmm^ 

Jo 2sinh(^) 
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l)|| 2 |+4|*| = p(l)|| 2 |( 



sinh 2 (|6>|/2) 



coth(|0|/2) 



+ 1 



Substituting the expression for the length, we get 



:i)in+4N 



2(\\9\ - sinh(|6»|)| + 2 sinh 2 (|#|/2)) ' 
In the Figure [3] we present the domain reachable from the origin by timelike geodesies. 



□ 





2.5 



Figure 3. The domain reachable from the origin by timelike geodesies 



In the case of spacelike curve we prove the analogous theorem 

Theorem 4. Let B = (x,y,z) be a point such that —x 2 + y 2 > 0, —1 < i < 1 ; then 

there is a unique spacelike geodesic, joining O = (0,0,0) with the point B. Let 9 satisfy the 
equation fj4.7j) Then the equations of spacelike geodesies are given by f|4.8j) and (|4.9p The square 
of the length of these geodesies is 

t 2 = 2 !M!±ilf! 

2(||fl| - sinh(|0|)| + 2sinh 2 (|^|/2)) ' 

Remark 2. Notice that (14. 4p and (|4.5p imply that the horizontal lightlike geodesies starting 
from the origin satisfy the equations x(t) = ±y(t), z(t) = 0. 

Remark 3. The result of Theorem [3] gives the description of the set reachable from the origin 

by timelike geodesies. Grochowski 111] proved that the region ^ a < 1, x > is reachable 

y 

by horizontal timelike future directed curves. Theorem [3] states that this set is reachable by 
geodesies of the same causal character. 
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5. HAMILTONIAN AND GEODESICS FOR H WITH THE Lorentzian METRIC 

In the present section we find the parametric equations of geodesies for quaternion H-type 
groups furnished with the Lorentzian metric. The horizontal wave operator is 

2A H = -X\ + X 2 2 + X\ + X\ 

Q 2 Q2 Q2 Q2 Q2 Q2 Q2 

X2X47; 7\ 1" X 2 X 3 - 7\ X 3 X 4 - 



dxf 0x2 ax§ ox\ OZ1OZ2 OZ1OZ3 OZ2OZ3 

1 d 2 1 d 2 1 d 2 

. 1 U / 2 _ 2 1 2 1 2\ , 1 V 1 2 1 2 1 2 _ 2\ , 1 _^ / 2,2 2 , 2\ 

~i~dz 2 %2 Xs X 4<9? X2 X3 X/l \~~dz 2 2 x 3 "r X 4J 

d ( d d d d \ d ( d d d d \ 

+ ~ X2 a h X 4 a X 3^— + x 4 - X 3 - h X27; h Xl-— 

azi \ axi 0x2 0x3 0x4 J 0Z2 \ ox\ 0x2 0x3 0x4 J 

d ( d d d d 

+ 7r~ X 3 o h X4T, V Xl- X2T. — 

OZ3 \ OXI OX2 0x3 0x4 
Then the associated Hamiltonian function H (£, 6, x, z) has the following form 

H = + &) + ^2X4^1^2 + X 2 X3e i 03 - X 3 X 4 # 2 # 3 ) 

+ ^9\{x\ -X 2 2 +X 2 3 + xl) + -6l{x\ + x\ + X% - xl) + -6l{x\ +X 2 2 -X 2 3 + x\) 

(5.1) +^i(-^26 - x l6 + x 4 £ 3 - £361) + 2^2(^1 - + x 2 & + xi£ 4 ) 

+ ^3(^36 + X 4 £ 2 + Xl6 - X 2 60- 

The corresponding Hamiltonian system is 
dH 

Xl = 7^- = -6 - 1x26*1 + ix 4 6* 2 + 5X3^3, 

X2 = TT— = 6 - 2: X 1^1 ~ l X 302 + 5X46*3, 
C?2 

dH 

X3 = TTT- = £3 + \x46l + tt22#2 + \x1O3, 
8^3 2 2 2 

dH 

X4 = tjt- = £4 - 5X36*1 + \x\d2 ~ \x263, 
dH 

dH 

Z2 = -QQ- = \{x 2 X49l ~ X3X463) + \0 2 {x\ +x\+x\~ xl) + \{X4^1 - X 3 £ 2 + X 2 £ 3 + Xl£ 4 ), 

dH 

^3 = 717- = \{x2X30i - X3X462) + l0 3 (xj +x\-x\ + xl) + i(x 3 £i + X4C2 + xi£ 3 - x 2 &t), 

Cf73 

6 = = -\xM + e 2 + e 2 ) + \^ei - \uo 2 - fae 3 , 

dH 

6 = ~d~X~ = "2 X2( "^ + 02 + °® + " ^° 2 + ^ ~ \ X ^° 2 ~ I X 3#1#3, 

dH 

^ = ~dx~3 = "2 X3( ^ + ° 2 ~ ^ + ^ + ^ ~ ^ 103 ~ \ X2 °^ + l X 4^2^3, 

£4 = = — I x 4(^i ~ ^2 + ^3) — i^3^i _ ^i^ 2 _ ^£2^3 - \x2Q\Q2 + ^Xa^^, 



ii = tttt- = 5(x 2 x 4 6'2 + X2X 3 6>3) + \0\{x\ -xl + x\ + xl) + \{-X2^1 - x^ 2 + x 4 £ 3 - x 3 £ 4 ), 
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02 = 



dH 

dz 2 
dH 

dz 3 



0, 
0, 

0. 



We observe that 61,62, 63 are constants. Let us remind that the projection of a solution of 
the Hamiltonian system onto (x, z)-space is called geodesic. In order to find it we will reduce 
the Hamiltonian system to the system containing only (a?i, x 2 , X3, £4, Z\, z%, Z3) coordinates. If 
we express £1, . . . ,£4 from the first 4 equations and substitute them in the equations of the 
Hamiltonian system, then we obtain 



6 



l^26i 



(-±i< 



X462 - X363), 



1 - X 3 6 2 + ±463) 



£3 = g-OMi + x 2 6 2 + xi6» 3 ), 

?4 = ^{-X2,6i + Xi6» 2 - X 2 6 3 ). 

Differentiating first 4 equations and substituting £1, ... ,£4 there, we get 

xi = -x 2 6i + X462 + £363, 

X.2 = -X161 - X362 + X463, 
X3 = X461 + X 2 6 2 + Xi& 3 , 
X4 = -X361 + X162 - X263 

or 



(5.2) 



We are looking for the solution x\ = x±(t), ... ,X4 = X4(t), t £ [—00, +00], satisfying xi(0) = 
0, . . . , #4(0) =0 and ii(0) = x\, . . . , £4(0) = £4. The eigenvalues of the matrix 

/ 







f 


— 6\ 63 62 






fx, \ 


X2 




-6l 


-62 63 






X2 


X3 




63 


62 0i 






X3 


\X4 J 




\ 62 


-03 -01 


J 




\X4 J 



A 







63 

V o 2 



-(71 



-0 3 







#2 \ 
63 
6l 

7 



are Ai = a, X2 = —a, A3 = ia, and A4 = — ia, where a = \J6\ + 6\ + 0|. 
eigenvectors are 

ui = (a, -01,03,02), 

v 2 = (a, 0i,-03, -02), 

^3 = (0, 0103 + ia6 2 , 6l + 6l, ia6 x - 2 3 ), 

V4 = (0, 0i0 3 - iad 2 , el + 0|, -ia0i - 2 3 ). 



The associated 
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The solution of the system (|5.2p is of the form 

x\ = a{c\ + c 2 ) cosh(ai) + a(c\ — 02) sinh(ai), 
%2 = 61(02 — ci) cosh(at) — 9\{c\ + 02) sinh(ai) 

+ 2(036163 — 04,0,62) cos(at) — 2(046163 + 03062) sin(ai), 
X3 = ^3(ci - c 2 ) cosh(at) + 63(01 + c 2 ) sinh(at) 

+ 2c 3 (6 2 + 6l) cos(at) - 2c 4 (6l + 6 2 2 ) sin(ai), 
X4 = 62(01 — c 2 ) cosh(at) + 6 2 (c\ + 02) sinh(ai) 

— 2(0406*1 + 036263) cos(ai) — 2(0306*1 — C46263) sin(ai), 

where 

ax? - 6 lX ° 2 + 63^3 + 6 2 ±l 

Cl = 2a~2 ' 

_ OX? + flixH - &3±l - e 2 ±l 

° 2 ~ 2o2 ' 

_ 6 X 6 Z ±\ + (gl + ^3 ~ 0293X1 

° 3 2a 2 (6 2 + 6 2 ) 

6 2 ±° 2 + flix^ 
C4 " 2a(6*2 + 61) " 



Therefore, 

xi = Ay sinh(at) + B\ cosh(at) + Ei, 

(5.3) X2 = A2 sinh(at) + B2 cosh(ai) + C2 sin(at) + D 2 cos(at) + E2, 
X3 = A3 sinh(a£) + B3 cosh(ai) + C3 sin(ai) + D3 cos(at) + E3, 
X4 = A 4 sinh(at) + B4 cosh(at) + C4 sin(at) + D4 cos(ai) + E 4 , 

where the coefficients Ai, Bi, Cj, Dj, Ei for i = 1,2, 3, 4, j = 1, 2, 3 are given in terms of q, 6j 
and are presented in Appendix. 

We find from the horizontality condition (|3.5p 

z\ = aa^t + Pi sinh(at) sin(at) + /3 2 cosh(ai) cos(ai) 
+/3 3 sinh(at) cos(at) + fi\ cosh(ai) sin(at) 
+a\ sinh(ai) + a\ cosh(ai) + a\ sin(at) + a\ cos (at) — 

(5.4) Z2 = aa^t + 01 sinh(at) sin(at) + 01 cosh(at) cos(at) 

+(3 2 sinh(at) cos(ai) + fi\ cosh(at) sin(at) 
+a\ sinh(ai) + ag cosh(ai) + a 2 sin(at) + ot\ cos (at) — 01, 
Z3 = aafjt + pf sinh(at) sin(at) + fi\ cosh(at) cos(ai) 

+P3 sinh(at) cos(at) + fi\ cosh(at) sin(at) 
+a| sinh(ai) + a\ cosh(ai) + Qy sin(at) + a| cos(ai) — (3 2 . 
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The coefficients and can also be found in Appendix. Then after long and laborious 
calculations we get 

(5.5) - x\ + x 2 + x\ + x\ = -16cic 2 sinh 2 (y) + 16£;sin 2 (y). 
and 

z \ + z 2 + z i = 4(at(— c\C2 + k) + 2c\C2 sinh(at) — 2k sin(ai)) 2 
-4k(cje at + c 2 e~ a *)(4sin(ai)sinh(at) + 5cos(at) - 5cosh(at)) 

(5.6) +8cic 2 &;(5 sin(ai) sinh(ai) +4cos(at) — 4cosh(ai)). 

Notice that — cic 2 + k = -^WvqW 2 . The auxiliary calculations can be found in Appendix. 
We believe that the quantity of geodesies in case of quaternion group equipped with the 

Lorentzian metric depends on the ratio , 3 1 3 2 z 3 ^ 2 like in the case of the Heisenberg group. 

\ x l x 2 x 3 X 4I 

The study of this question we postpone for a forthcoming paper. 



6. Appendix 
The coefficients for the solutions (15.31). 



A l = c l + C2 , A 2 = d -^p^, A, = °^-^, M= 9 -^p±. 

o o gl(ci+C 2 ) ^(^+€2) ^(cA^) 

-Dl = Cl — C 2 , i3 2 = , -D3 = ; £>4 = • 

a a a 

2(c 3 0i0 3 - c 4 a# 2 ) ~ 2c 3 (# 2 + # 2 ) „ 2Mi + c 3 M 3 ) 

t>2 — j ^3 — , W 



n 2(04^1^+ c 3 a0 2 ) 2c 4 (#?+# 2 ) 2(c3O0i -^Ms) 







„ „ Oi(ci+c 2 ) 2(c 4 0i03 + c 3 a0 2 ) 

E 1 = c 2 -ci, E 2 = - 



# 3 ( Cl +c 2 ) 2c 4 (0 2 + 0|) _ 2 (ci + c 2 ) 2(c 3 4-qM 3 ) 

£/ 3 — , £,4 — . 

a a a a 

The coefficients and in © are related in the following way. f3\ = + a\), 

P2 = -r{— oc\ +a\), (5\ = -(— a 2 + a\), (3\ = -(a 2 + a 3 ), i = 1, 2, 3. The list of is presented. 

«o = ^(-cic 2 + fc), 
a 2 = 2 &(_ ClC2 + A: ) ) 

«o = ^r(-ciC2 + k), 

a \ = |( c i + c 2)( c 3^i6 | 3 - ac 4 6> 2 ), 

Q i = |(ci + c 2 )(c 3 6i 2 6i 3 + ac 4 6»i), 

«? = -^(ci + c 2 )(^ 2 + ^), 
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„1 _ 4, 



a\ = -(ci - c 2 )(c 3 0i0 3 - ac A 9 2 ) 

a \ = a( c i ~ c 2 )(c 3 6 2 6 3 + ac 4 0i) 

«1 = -^(ci-c 2 )(0f + 2 2 ), 

a 3 = |( c i + C2)(c46'i6 | 3 + ac 3 6> 2 ) 

Q 3 = a( C l + C 2)(C46 , 203 ~ ac 3 0i) 



«S = -^(ci + C2 )(^+0 2 2 ), 

«4 = | ( c l ~ C 2 )(c 4 0i0 3 + 0036*2), 

«4 = l(ci - c 2 )(c 4 2 03 - ac 3 0i), 
«1 = -^(ci-c 2 )(0 2 + 2 2 ), 

a* = -20 2 (c 3 ( Cl + c 2 ) - c 4 ( Cl - c 2 )) + ^Siafi - ^(c 3 ( Cl - c 2 ) + c 4 ( Cl + c 2 )), 
a 2 = 2^(03(0! + c 2 ) - c 4 ( Cl - c 2 )) + ^1 _ «3 (C3 ( C1 _ C2 ) + C4 ( Cl + C2)); 
«l = + # 2 2 )(c 3 (ci - c 2 ) + C 4 ( Cl + c 2 )) + 

a\ = 20 2 (-c 3 (ci - c 2 ) + c 4 (ci + c 2 )) - ^(c 3 (c! + c 2 ) + c 4 (ci - c 2 )), 
a§ = 20i(c 3 (ci - c 2 ) - c 4 (ci + 02)) - ^(c 3 (ci + c 2 ) + c 4 (ci - c 2 )), 
«6 = s(^i + ^)(c 3 (ci + 02) + c 4 (ci - c 2 )), 

a l = -20 2 (c 3 ( Cl + c 2 ) + c 4 ( Cl - c 2 )) - 2^(-c 3 (ci - c 2 ) + c 4 (ci + c 2 )) - ^ife, 
a 2 = 20 1 (c 3 (c 1 + c 2 ) + c 4 (ci - 02)) - ^(-c 3 (ci - c 2 ) + c 4 (ci + c 2 )) - 1 2 fc, 

«7 = ~Wl + ^2 2 )(C 3 (C1 - C 2 ) - C 4 ( Cl + C 2 )) - ±0 3 k, 



1 12 2 3 3 

a 8 = _a 6' a 8 = _a 6' a 8 = _a 6' 



where k = (c| + cl)(0 2 + 0|). 

In the following lines we give some useful calculations to deduce (|5.5p . 



-A\ + A 2 + A 2 + A 2 = -4ciC2, --B 2 + £ 2 + £ 2 + £ 2 = 4cic 2 , 
C 2 + C| + C| = D\ + D\ + D\ = 4k, -E\ + £ 2 + £| + E 2 = 4{ Cl c 2 + k) , 

-AxBx + A 2 B 2 + A 3 B 3 + A A B A = A 2 C 2 + A 3 C 3 + A A C A = A 2 D 2 + A 3 D 3 + A A D A = 0, 

-A X E X + A 2 E 2 + A 3 E 3 + A A E A = 0, 
B 2 C 2 + B 3 C 3 + B A C A = B 2 D 2 + B 3 D 3 + B A D A = 0, 
-B X E X + B 2 E 2 + B 3 E 3 + B A E A = -4cic 2 , 
C 2 D 2 + C 3 D 3 + C7 4 L> 4 = C 2 £ 3 + C 3 E 3 + C A E A = 0, D 2 E 2 + L> 3 £ 3 + D A E A = -4k. 
Now we present the values of different combinations that we used in calculation of z\ + z\ + z 2 . 

(<*l? + (a 2 ? + (a 3 ? = 4(-c lC2 + k) 2 = ^(-(4? + (x°) 2 + (x°) 2 + (±°) 2 ) 2 , 

(A 1 ) 2 + (Pi) 2 + (Pi? = (Pi) 2 + (Pi? + (/3 2 3 ) 2 = (Pi? + (Pi? + (Pi? 

= (p 1 A ? + (P 2 A ? + (P 3 A ? = 2(c 2 + c 2 )k, 

(alf + (a 2 ) 2 + (al) 2 = 8(c? + c 2 )k + 16c 2 c 2 , 
(ai) 2 + (a 2 ) 2 + (a|) 2 = 8(c 2 + C 2 )A ; , 
(a*) 2 + (a 2 ) 2 + (a?) 2 = 8(c 2 + c 2 )fc + 16fc 2 , 
(al) 2 + (a 2 ) 2 + (al) 2 = 8(c? + c 2 )k, 

alp\ + a\0[ + = a 1 ^ + a 2 (3 2 + a 3 ^ = a^ 1 + a 2 ,/? 2 + c^f = a^ 1 + ag/?| + ag^ = 0, 

ajag + «o q2 + a o a l = — ^c 2 c 2 + 8cic 2 /c, a o a 6 + a o a 6 + a o Q 6 = 
a^aj + «q0 2 + aga^ = 8cic 2 k — 8k 2 , a o a l + a o a i + a o a l = ^> 
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PlPl + p\p\ + + PlPl + p!p! + PlPi = o, 
PlPl + P 2 Pi + pfPl = o, 
PlPl + PlPl + p!p! = 2(c? - 4)k, 

alPl + a 2 /? 2 + alPf = -4(c 2 - c 2 2 )k, 
a\0[ + a\(3\ + + a 2 /? 2 + ^ + ^3 = _ 8(c 2 + C | )A . 

a i/?i + a 2^ + a^s = o, a^g 1 + a\(3\ + a 2 -/? 2 + a 2 /? 2 + a?/3g + a |/3? = 0, 

(31(31 +(3 2 2 (3l + PlPl = 2(cl-c 2 2 )k, (3^l+PlPi + PiPl = 0, 

a\(3\ + a\(3\ + a\i3 2 + a 2 (3 2 + «i/3 2 3 + a^f = 0, 
a\fi\ + ag/3 2 + a lPl = 8cic 2 fc, 
+ a\(3\ + a 2 7 (3l + a|/9| + c? 7 p\ + af/3f = 0, 

ot\(3\ + a\(i\ + a\(3\ = -8c\c 2 k, a\(3\ + a 2 (3 2 + ag^f = -8cic 2 /c, 
al(3l + a 2 P 2 + a|/?f = 0, a\{3\ + a 2 /? 2 + a\l3\ = -4(c 2 - c 2 )k, 

a\p\ + a?/?| + a 3 7 (3l = -8ac 2 k, 

a\a\ + a\ a\ + OgOg = 8(cf - c^fc, a^Oy + a\a 2 + afay = 0, 
a\a\ + a§ + a|a| = -8(c 2 - c 2 )fc, aga 7 + a\a 2 + agaf = 0, 
a\a\ + aga| + a%a\ = — 8(cf + c 2 )/c, a\a\ + a 2 ctg + a 7 ag = 0, 

P\Pl + /3 2 /3 2 + plPl = -4 Cl c 2 k, 
alPl+alPl+ulPl = 0, ^Pl + c^pl + c^pl = -4(c 2 -c 2 )/c, a\p l 2 +u 2 p 2 +al(3l = -8 Cl c 2 k. 
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